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The peculiar features of quantum magnetism sometimes forbid the existence of gapped ‘featureless’ 
paramagnets which are fully symmetric and unfractionalized. The Lieb-Schultz-Mattis theorem is 
an example of such a constraint, but it is not known what the most general restriction might be. 

We focus on the existence of featureless paramagnets on the spin-1 square lattice and the spin-1 
and spin-1/2 honeycomb lattice with spin rotation and space group symmetries in 2-l-lD. Although 
featureless paramagnet phases are not ruled out by any existing theorem, field theoretic arguments 
disfavor their existence. Nevertheless, by generalizing the construction of Affleck, Kennedy, Lieb 
and Tasaki to a class we call ‘slave-spin’ states, we propose featureless wave functions for these 
models. The featureless-ness of the spin-1 slave-spin states on the square and honeycomb lattice are 
verified both analytically and numerically, but the status of the spin-1/2 honeycomb state remains 
unclear. 


I. INTRODUCTION 

In the Landau paradigm phases of matter are classi¬ 
fied by their symmetries, and fully symmetric phases are 
ubiquitous. Only after spontaneous breaking of these 
symmetries do classical orders emerge. One reason the 
Landau paradigm is incomplete is that a gapped and fully 
symmetric phase is sometime excluded in a zero temper¬ 
ature quantum phase diagram. For example, the Lieb- 
Schultz-Mattis (LSM) theoreirfi] asserts that an S=l/2 
spin chain in 1-l-ID does not admit a gapped paramag¬ 
netic phase that respects all the symmetries, specifically 
translation and 50(3) spin rotations; the system is either 
gapless or breaks a symmetry. A generalized version of 
LSM theorem in higher dimensions finds that a magnet 
with half-integer spin per unit cell cannot be a short- 
range entangled paramagnet; the system is either gap¬ 
less, breaks a sym metr y, or has fractionalized excitations 
(topological order) .1^ We refer to a gapped, symmetric 
and short-range entangled paramagnet as a ‘featureless’ 
paramagnet. The LSM theorem only requires transla¬ 
tion and spin rotational symmetries, though versions in¬ 
corporating additional space-group symmetries are also 
possible.03 

To determine whether a featureless phase can exist in 
a model (as defined by its Hilbert space and symmetries 
- not a Hamiltonian), one must either prove an extension 
of the LSM theorem forbidding its existence, or produce 
an explicit example. Here we focus on the spin-1 square 
lattice and the spin -1 and spin- 1/2 honeycomb lattices in 
2-l-lD. Although not forbidden by any existing theorems, 
it is unknown whether featureless phases exist for these 
models. In fact, field theoretic arguments based on the 
non-linear sigma model find non-trivial spatially depen¬ 
dent Berry phase pattern^SHU which strongly disfavor the 
existence of such featureless phases, at least proximate 
to Neel or valence-bond orders where the field theory is 
justified. 


Finding featureless paramagnetic wave functions is 
non-trivial as they require a finite amount of quantum 
entanglement. For example, Ref. IH] first explored this is¬ 
sue for the spin- 1/2 honeycomb model at zero magnetiza¬ 
tion, requiring only U(l) symmetry, and found that only 
strongly interacting states can provide an example of a 
featureless paramagnet Requiring full SO (3) symmetry 
further complicates the search. In this work, we propose 
and analyze a class of featureless SO(3)-symmetric wave 
functions for the square and honeycomb magnets, and 
find conclusive evidence for the existence of featureless 
spin -1 paramagnets on the square and honeycomb lat¬ 
tices. Our numerics are unable to establish whether the 
proposed spin- 1/2 honeycomb state has exponentially de¬ 
caying correlations. 

The wave functions we propose are written in terms 
of projected entangled pair states (PEPS^P^^. In gen¬ 
eral, global symmetries impose strong constraints on the 
PEpglMill, To serve our purpose, we choose the ‘en¬ 
tangled pairs’ to be multiplets of SO(3) symmetric spins 
and, thus, will refer to them as the ‘slave-spin’ PEPS. 
By construction these states are fully symmetric, though 
numerical evidence is required to prove that they are 
gapped states without spontaneously broken symmetry 
(“cat states”) or topological order. These slave-spin 
PEPS suggest a generalized version of the deconfined 
quantum criticality when brought to a second order phase 
transition to a classically-ordered phase. 

The rest of the paper is organized as follows. In Sec. |H] 
we discuss the ID spin-1 chain as a warm up, in order to 
demonstrate in a simple setting how field theoretic argu¬ 
ments based on non-trivial Berry phases (the topological 
0-term) can be evaded through the slave-spin PEPS con¬ 
struction. In Sec. m we study the spin -1 square lattice. 
After reviewing the field theoretic obstacl^^ to a fea¬ 
tureless state, we construct a featureless wave function 
with spin rotation, translation, C 4 lattice rotations, and 
time reversal symmetry. A numerical analysis verifies its 
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featureless-ness. In Sec. |IV| we study the spin-1 and 
spin-1/2 honeycomb lattices. We derive a similar field 
theoretic obstacle for these models, and construct a wave 
function with SO(3) spin rotation, translation, C 3 lat¬ 
tice rotations, mirror, and time reversal symmetry. For 
the spin-1 magnet, our numerics unambiguously verify its 
featureless-ness, while for spin-1/2 it is unclear whether 
the proposed state is gapless, or simply has a large corre¬ 
lation length we cannot yet detect. In Sec. |V]we conclude 
with a discussion of open questions and future directions. 


II. 1-blD FEATURELESS PARAMAGNETS ON 
THE SPIN-1 CHAIN 


We warm up with a 1-l-lD spin chain with spin rotation 
and translation symmetry, both to review the field the¬ 
oretic Berry-phase analysis and introduce the slave-spin 
formalism in a simple setting. The field theory proceeds 
by deriving an effective action for the Neel order param¬ 
eter n S .1^ For a spin chain with spin S per site, in 
addition to the kinetic terms there is a non-trivial topo¬ 
logical Berry phase (the ©-term^PS; 

0 f 

S'top[n] = ^ / dxdrn ■ {dj;n x dru), (1) 

where 0 = 27rS. For integer spin there is a gap in the 
bulk, leading to a featureless gapped paramagnet, the 
Haldane phase.For 5 = 1 the term 5top[n] nevertheless 
leads to interesting degenerate edge states for an open 
chain. At S = 1, 5top is a total derivative in the bulk, 
but induces a 0(3) Wess-Zumino-Witten (WZW) term 
at level 1 on the boundary, which implies a spin-1/2 de¬ 
gree of freedom is localized on the boundary of the spin 
chaiiP. 

An ideal wave function for the spin-1 Haldane phase 
and its edge states is given by the Afheck-Kennedy-Lieb- 
Tasaki (AKLT) state, ^ the simplest example of a ‘slave- 
spin’ state. We start by associating two auxiliary spin- 
1/2 degrees of freedom labeled by \ai(3i) on each site 
i, with ai,f3i =t / /. A featureless state |'I'o) for the 
auxiliary spins is obtained by projecting the |/3i) on the 
Ah site and |ai+i) on the (i -\- l)th site onto the spin-0 
channel: 


1 ^ 0 )= E n I I ...), (2) 

{a,P} \ i ) 


with 5 a 2 X 2 matrix: 

5 = 


0 1 
-10/' 


(3) 


S is obtained from the Clebsh-Gordon (CG) coefficient 
of “fusing” two spin-1/2s into a spin singlet, namely 
Safi = (S = 0|S = 1/2, a; S = l/2,/3). To obtain a 
wave function for the physical spin-1 degrees of freedom, 
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Figure 1: (a) represents the projector of two spm-1/2 

degrees of freedom onto the spin singlet channel. The arrow 
indicates the anti-symmetric property of Sap- (b) Tj)/ with 
symmetric subscripts represents the projector of two spin- 
1/2 degrees of freedom onto the spin triplet channel, (c) The 
graph represents the tensor network the AKLT wave function. 
The red lines represent physical spin-1 degrees of freedom, 
whereas the joined black lines represent the auxiliary spin- 
1/2 degrees of freedom that are contracted. 


(b) 
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we then project the two auxiliary spin-l/2s onto a spin-1 
channel of each site: 


^AKLT({mi}) = (...mjmi+i...|4')o 

= E ’ (4) 

{a./3} V * / 


where mi = 1,2,3 represent the three spin-1 states on 
the Ah site. The tensor = (S = l,m|S = 1/2, a; S = 
1/2, fi) is the GG coefficient for fusing two spin-l/2s into 
a spin triplet. Graphical representations of the tensors 
Sap and and their contraction are shown in Fig. [T 
The wave function of Eq. is actually in the form o ’ 
a matrix product state (MPS), as it takes the form of 
a matrix-trace in the auxiliary a/(3 space. The ‘bond 
dimension’ x of the MPS is y = 2, corresponding to two 
the auxiliary spin states. For a finite open spin chain 
each boundary has one un-summed matrix index that 
corresponds to the spin-1/2 edge state of the Haldane 
phase, consistent with the field theory prediction. 

Is the boundary spin-1/2 inevitable in any paramag¬ 
netic phase of the spin-1 anti-ferromagnet? No: it is sim¬ 
ple to construct a slave-spin counterexample to this field 
theoretic expectation. Instead of the auxiliary spin-l/2s 
of the AKLT case, we start with two auxiliary spin-ls on 
each site f, with the states \lini),li,ni S 1,2,3 labeling 
the Sz basis. We again project pairs \ni) and \li+i) (for 
all i) into the spin-0 channel by using the GG coefficient 
Snji^i = (S = 0|S = l,ni;S = l,^i+i). We then project 
the auxilliary spins of each site Ihrii) onto the total spin-1 
channel in order to obtain a spin-1 chain wave function. 
This projection is accomplished by the GG coefficients 
= (S = l,m|S = l,^i;S = l,ni). The new spin-1 
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Figure 2: (a) Sap is a symmetric matrix that represents 

the projector of two spin-1 degrees of freedom onto the total- 
spin-singlet channel, (b) with anti-symmetric subscripts 
represents the projector of two spin-1 degrees onto the spin- 
1 channel, (c) The graph represents tensor network for the 
wave function l^'para)- The red lines represent physical spin-1 
degrees of freedom, whereas the joined black lines represent 
the auxiliary spin-1 degrees of freedom that are contracted. 


under the following symmetries: the global S'0(3) spin 
rotation, the translation, the lattice C 4 rotation (i?^) 
the time reversal (T) and the mirror symmetries (dv with 
the vertical mirror plane and 17450 with the 45° mirror 
plane). We will first review the field theoretic obstacle to 
such a spin-1 featureless paramagnetic phase. Then, we 
prove its existence by providing an example wave func¬ 
tion using the slave-spin PEPS construction. A numer¬ 
ical analysis verifies its featureless-ness. In the end, we 
will comment on the connection to previous studies Ref. 
\TmE\ on the spin -1 square lattice. 


slave-spin wave function is 

l^-para) = I • • • (5) 

{m} \ i / 

The graphical representation of the tensors S and T is 
given in Fig. Since the state is composed by suc¬ 
cessive SO(3) invariant projections, by construction it is 
invariant under global spin rotations and is translation 
invariant. However, we must verify it is not a ‘cat-state’ 
of symmetry breaking valence-bond patterns. The stan¬ 
dard MPS test for such states is to analyze the ‘transfer 
matrix’ (Fig. (d)) defined by 

= ( 6 ) 

m 

If the dominant eigenvector of Inn',u' is unique, which we 
find it is, all equal-time correlation decay exponentially, 
so I'Fpara) IS a featureless paramagnet. 

Since we used spin -1 auxiliary degrees of freedom, on 
an open chain the edge state of iTpara) are effectively 
spin-1, counter to the field theoretic expectation. In fact, 
l^para) represents the “trivial phase” in the Z 2 classi¬ 
fication of SO{3) spin symmetry protected topological 
states, while the AKLT state represents the “non-trivial 
phase”l^3H2ll, 

In conclusion, while the seemingly universal field theo¬ 
retic analysis is highly suggestive (the ground state of the 
nearest-neighbor Heisenberg chain is indeed in the Hal¬ 
dane phase) it does not ‘mandate’ any properties of the 
energy spectrum in the same sense that the Lieb-Schultz- 
Mattis theorem does. 


III. 2+lD SPIN-1 FEATURELESS 
PARAMAGNET ON THE SQUARE LATTICE 

In this section, we study the featureless paramagnetic 
phase on the square lattice with spin 1 per site in 2-l-lD 


A. Field Theoretic Obstacle to Featureless 
Paramagnetic Phases on the Square Lattice 


From a field theory point of view, to obtain a param¬ 
agnetic phase, one can start from an classically ordered 
background and disorder it through the proliferation of 
topological defects/excitations or instantons. In the Neel 
background in 2-l-lD, the topological instanton to be con¬ 
sidered is the monopole tunneling eveniP which happens 
at plaquette centers and changes the winding number 
(or the homotopy class) of the spatial configuration of 
the Neel order parameter at a given time (see Fig. 

(a) and (b)). There are non-trivial location-dependent 
Berry phase associated to the monopole tunneling events 
in the spin-coherent-state path integraP^. With spin S 
per site, this Berry phase is (-1-1)^®, (—1)^® or (—i)^® 

depending on whether the dual lattice coordinate (cc, y) 
of plaquette center where the monopole tunneling event 
happens is (even,even), (even, odd), (odd, odd) or (odd, 
even) (see Fig. 1 (c)). 

For S = 1/2, this non-trivial Berry phase leads to the 
deconfined quantum criticalitj®. When the monopole 
tunneling event is proliferated in the spacetime and the 
Neel background is disordered, the non-trivial location 
dependence of this Berry phase in Fig. |^(c) leaves an im¬ 
print, either breaking a space group symmetry or result¬ 
ing in topological order. Therefore, a featureless param¬ 
agnetic phase does not exist for S = 1/2 on the square 
lattice, which agrees with the generalized LSM theorerrP. 
Similar results applies for S = 3/2. For S = 1, even 
though the generalized LSM theorem does not apply, the 
Berry phase still has a non-trivial spatial pattern which 
results in either classical orders, a nematic phas^ for 
instance, or topological state^^ when monopole tunnel¬ 
ing events are proliferated. The obstacle to a featureless 
paramagnetic phase remains for S = 1. Furthermore, the 
Berry phase pattern is trivial when S = 2, in which case 
an example of featureless paramagnets is given by the 2D 
AKLT states on the square lattice. 
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Figure 3: A spatially uniform configuration of the Neel order 
parameter with a trivial winding number (a) tunnels into a 
configuration in a different topological sector with winding 
number 1 after a monopole tunneling event occurred at the 
plaquette center, (c) shows the spatial dependence of the 
non-trivial Berry phase associated to the monopole tunneling 
event in a spin-S Neel background. 


B. Spin-1 Featureless Paramagnets on the Square 
Lattice 

Does the field theoretic obstacle strictly forbid the ex¬ 
istence of a spin-1 featureless paramagnet? The answer is 
no: we can construct a slave-spin PEPS which serves as 
an example of a spin-1 featureless paramagnet. We start 
with 4 auxiliary spin-1/2 degrees of freedom per site on 
the square lattice with each spin 1/2 attached to a dif¬ 
ferent link connected to the same site. On each link, we 
project the two spin l/2s into a singlet pair using Sa/s- 
The physical spin-1 function Idfi)) is obtained by a second 
projection that project the 4 auxiliary spin 1/2 degrees 
of freedom on each site to a total-spin-1 channel. In gen¬ 
eral, the projection from 4 spin-1/2 degrees of freedom 
to a total-spin-1 channel is not unique. However, there is 
only one choice that guarantees the , T, Uv and CT45CJ 
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symmetries. Its graphical representation is shown in Fig. 

where we first form triplet pairs between the upper 
and lower spin l/2s and between the left and right spin 
l/2s and project the two triplet pairs into a total-spin-1 
channel. The associated tensor is given by 

jjm _ \ ' rpm rjinirpn-i /y'l 

^nesw / ^ ^ nin2 ^ ew ^ ns ' \ ) 

ni,2 

where ni ^2 sum over 1, 2 and 3 (which label the interme¬ 
diate spin-1 states) and the tensor T and T are the CG 
coefficients previously defined. The graphical represen¬ 
tation of I'I'u) is now given by Fig. 

I'!'£)) is by construction translationally invariant. Since 
all of the tensors involved are the CG coefficients, the 
global SO(3) spin rotation symmetry is also preserved. 
Under the site-centered or plaquette-centered (say 

counter-clockwise), each tensor contributes a — 1 

phase factor because under the cyclic permutation of the 
subscripts 

jjm _ _ jjm /o'! 

On top of that, the arrow on every horizontal link is 
flipped after the 7r/2 counter-clockwise rotation, each 
producing a —1 factor, while all the vertical ones are 
intact. Notice that the square lattice in a planar or torus 



Figure 4: The graphical representation of the tensor 
that implements the projection of 4 spin-1/2 degrees of free¬ 
dom to a total-spin-1 channel compatible with Rzjl and T. 
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Figure 5: (a) This graph represents the tensor network of the 
slave-spin PEPS Itfo). There is a singlet projector Sap on ev¬ 
ery link and tensor on every site, (b) shows the graph¬ 

ical representation of the horizontal transfer matrix (HTM). 


geometry has equal number of sites and horizontal links. 
The state I'l'^)) is invariant under R 2 jr (with 0 lattice an- 
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gular momentum). Since both and Sap are real, 

|4'£i) is T-invariant. We can also verify that pre¬ 

serves tTv and cr 45 o. 

At this point, we have shown that I'I'u) preserved all 
the required symmetries as a PEPS wave function. Now 
we verify numerically that (1) no symmetry is sponta¬ 
neously broken (“cat state”) and (2) no topological order 
is present in in I'l'^)). We consider I'I'd) on a cylindri¬ 
cal geometry with finite circumference L and view it as 
an MPS. As long as the transfer matrix of this MPS 
has a unique most dominant eigenvector (MDE) and the 
physical correlation length stays finite as L —>■ 00 , we 
can ensure the absence of spontaneously broken symme¬ 
tries. Two types of cylindrical geometries with either the 
y direction or the 45° direction compact should be con¬ 
sidered, which leads to two types of transfer matrices: 
the horizontal transfer matrix (HTM) and the diagonal 
transfer matrix (DTM). We perform the exact diagonal- 
ization to the HTM and DTM for L = 2,3, ...,12. In 
each case, we obtain a unique MDE. The physical corre¬ 
lation length ^ is plotted against 1/L in the left panel of 
Fig. [^and found to saturate as 1/L —)■ 0. The unique¬ 
ness of the MDE for both transfer matrices and the fi¬ 
nite correlation length ensure the absence of the sponta¬ 
neous symmetry breaking for the ones compatible with 
either of the transfer matrices, which are the transla¬ 
tion, the SO{3) spin rotation, T, CTv and (7450. To verify 
the i? 22 r symmetry, we calculated the nearest neighbor 
(NN) spin-spin correlation function along the 4 differ¬ 
ent directions. They agree with each other at the value 
{S-S)^^ = —0.462 even though the cylindrical geometry 
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Figure 6: |$aklt) (I'f’AKLT)) i® obtained by stacking ID 
AKLT chain along the y (a:) direction. Both of them are 
nematic paramagnets without topological order. The slave- 
spin PEPS I^'d) can be obtained by stacking the two states, 
I^aklt) 8-iid I'f’AKLx)" each other and project the 

two spin Is per site into a total-spin-1 channel. 


Figure 7: Left panel: The physical correlation length ^ in the 
unit of lattice spacing is plotted against 1/L. The result is 
obtained from the DTM with L = 4, 6,8,10,12. The FITM 
result is similar. Right panel: Entanglement entropy of the 
semi-finite cylinder is plotted against the circumference L. A 
linear fit captures the “area law” and the TEE. 


is incompatible with Rzjl ■ Therefore, we expect Rzjl is 

4 4 

not spontaneously broken in the infinite-plane limit. All 
the symmetry properties of |^'d) are now confirmed. 

Moreover, the R^^ symmetry of \'^d) and the unique¬ 
ness of the MDE for both HTM and DTM also implies 
the absence of topological orders. The argument is the 
following. Let’s first assume that Id/o) is topologically 
ordered. We can construct the wave function It'd) on a 
X Ly torus with periodic boundary conditions along 
along the x and y direction. When = Ly, the R'^ 
symmetry of Iffr)) implies that |d>£)) is an eigenstate of 
the topological S matrijPSI An eigenstate of the topo¬ 
logical S matrix does not have a fixed anyon flux along 
either the x oi y cycle of the torus. In other words, 
the state on the torus is not a minimally entan¬ 

gled state (MESj^ along either directions. Therefore, 
the Wilson loops around the y {x) direction cycle de¬ 
velops a long-range order when we view the system as 
a quasi-ID system in the x (y) direction. If the system 
is gapped, this long-range order will persist even when 
we open up the torus to create a infinite cylinder (with 
fixed circumference). The HTM will develop degener¬ 
acy in its MDEs. Similar argument can be made for the 
DTM. However, this is inconsistent with the numerical 
result which shows a unique MDE. Therefore, we prove 
by contradiction that the PEPS Idfri) is not topologically 
ordered. To further justify this argument, we calculate 
the entanglement entropy (EE) of a semi-infinite cylin¬ 
der and plot it against the circumference L in Fig. We 
obtain an “a rea la w” plus a topological entanglement en¬ 
tropy (TEE) 21 H 2 | 7 ~ 0 , which indicates the absence of 
topological order. Therefore, we conclude |4'£)) is indeed 
a featureless paramagnet. 

Now, we have confirmed that |d>£)) is a featureless 
paramagnetic state which lies beyond the field theoretic 
expectation. Consequently, if there is a second-order 
transition separating the featureless paramagnet 
and the Neel order or the VBS order where the Berry 
phase pattern in Fig. is obtained, the critical the¬ 
ory should be a generalized version of the deconfined 
quantum criticality and most likely one that properly 
captures the “slaved” degrees of freedom in the slave- 


spin PEPS. The construction of such a critical theory is 
beyond the scope of this paper and will be left for fu¬ 
ture works. Apart from the field theory, one can also 
try to search for simple Hamiltonians such that 14 'd) is 
an approximate ground state. In the numerics, we cal¬ 
culate the nearest neighbor (NN) and the next nearest 
neighbor (NNN) spin-spin correlation functions and ob¬ 
tain {S ■ <S)nn = —0.462 and {S ■ 5)nnn = —0.174 Since 
both of them are anti-ferromagnetic (AF), a Hamiltonian 
with AF Ji and J 2 Heisenberg couplings should be prefer¬ 
able. The actual J 1 -J 2 Heisenberg model at intermediate 
J 2 /J 1 is shown to host a nematic phas^^that is captured 
by stacked AKLT chain^^. We denote the 2 symme¬ 
try related stacked-AKLT-chain states as |d>AKLT) 

I ‘^akltL Their tensor network representations are given 
in Fig.^ When we stack |$aklt) top of Id^AKUT) 
project the two spin 1 degrees of freedom on each site in 
to a total-spin -1 channel, we can obtain the featureless 
paramagnet I'I'/j) (see Fig. [ 6 ). It is an interesting ques¬ 
tion whether this connection between the wave functions 
leads to a dynamical theory that captures both the ne¬ 
matic phase and the featureless phase. 


IV. 2-l-lD FEATURELESS PARAMAGNETS ON 
THE HONEYCOMB LATTICE 

In this section, we will study the featureless param¬ 
agnets for spin- 1/2 and spin -1 honeycomb lattice with 
the following symmetries: the global spin rotation, the 
translation, the lattice C 3 rotation (R^), the time rever¬ 
sal (T) and the 2 link-centered spatial mirror reflections 
symmetries {cr± with the mirror plan perpendicular and 
cry parallel to the link). We will first review the field the¬ 
oretic obstacle to such featureless paramagnetic phases. 
Then, we prove the existence of the spin-1 featureless 
paramagnet by providing an example wave function in 
the slave-spin PEPS form. A numerical analysis verifies 
its featureless-ness. For the spin-1/2, we proposed two 
types of slave-spin PEPS with bond dimension y = 2,6. 
The PEPSs with x = 2 are found to be “almost feature¬ 
less paramagnets” with time reversal symmetry and at 
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Figure 8: This plot shows the spatial dependence of the 
non-trivial monopole tunneling Berry phase in a spin-S Neel 
background on a honeycomb lattice. 



Figure 9: The figure is the graphical representation of Eq. 

The directly connected lines represent 5-functions between 
the indices on their end points, while the lines decorated by a 
black dot represent the singlet projector S defined previously. 
The physical degrees of freedom are represented by the red 
lines and auxiliary ones by the black lines. 


least one of the mirror reflection symmetries broken. Last 
but not least, we proposed the x = 6 slave-spin PEPS for 
spin-1/2 honeycomb system. We show analytically that 
the slave-spin PEPS is fully symmetric. The numerical 
test for a gap is inconclusive, so we can only conjecture 
that the proposed x = 6 slave-spin PEPS is a featureless 
spin-1/2 paramagnet on the honeycomb lattice. 


A. Field Theoretic Obstacle to Paramagnetic 
Phases on the Honeycomb Lattice 


Similar to strategy described in Sec. Ill A we start 
with the classical Neel background with spin S per site 
and consider the proliferation of monopole tunneling 
events in the spacetime to obtain a paramagnetic phase. 
A generalized spi n-coh erent-state path integral or dimer 
model calculatiorP^^I23| shows that the monopole tunnel¬ 
ing events on the honeycomb lattice have a non-trivial 
location dependent Berry phase pattern shown in Fig. 
When the Neel background is disordered, this non-trivial 
Berry phase pattern rules out the possibility of spin-1 
and spin-1/2 featureless paramagnets. In contrast, when 
S = |Z, the Berry phase pattern becomes trivial. Indeed, 
the most natural the AKLT state (a featureless param¬ 
agnet) on the honeycomb lattice requires S = 3/2^. 


B. Spin-1 Featureless Paramagnets on the 
Honeycomb Lattice 

We prove the the existence of the spin-1 featureless 
paramagnet by writing down an example in the slave- 
spin PEPS form. We start with 3 auxiliary spin-1 de¬ 
grees of freedom on each site with each spin-1 degree of 
freedom associated to each link that connects to the same 
site. Then we project the two spin-1 degrees of freedom 
associated to the same link into a total-spin-0 channel. 
This projection is implemented by the tensor Sni defined 
previously. Then, the three auxiliary spin-1 degrees of 
freedom on each site are projected onto a total-spin-1 
channel to obtain a physical spin-1 wave function. The 
most symmetric choice of projection is given by the ten¬ 


sor 


~ ^abSmc + Sac^mb + Sbc^n 


(9) 


whose graphical representation is given in Fig. Here, 
the indices a, b, c and m all run from 1 to 3 each rep¬ 
resenting 3 spin-1 states in the basis. The directly 
connected lines represent S functions between the indices 
on their end points, while the lines decorated by a black 
dot represent the singlet projector S defined previously. 
This tensor is invariant under the permutation of all the 
subscripts and T: 


TTim* _ _ ZT’’^ _ Z?’^ 

^ahc ^abc -^bca ^acb’ 


( 10 ) 


After the second projection, we obtain the featureless 
paramagnetic wave function whose tensor network 

is shown in Fig. 10 (a). This tensor network has the 


tensor Sni on each link and the tensor on each site. 

Because the tensors and Sab are the CG coeffi¬ 
cients that are real and invariant under the permutation 
of their subscripts, the wave function I'kf) is invariant 
under all symmetries required. Now, we present the nu¬ 
merical evidence that I'kf) is free of spontaneously bro¬ 
ken symmetry and topological order. Due to its long 
correlation length, we cannot perform exact diagonaliza- 
tion study on finite-size cylinder. Instead, we directly 
study the planar limit. We view the wave-function norm 
('I'f I'I'f) shown in Fig. 10 (b) as the partition function 
of a classical statistical model for the auxiliary spin-1 de¬ 
grees of freedom. Any spontaneous symmetry breaking 
manifests itself as the non-uniqueness of the MDE for its 
transfer matrices. For the honeycomb lattice, there are 
two types of transfer matrices: the armchair transfer ma¬ 
trix (ATM) (Fig. [^(c)) and the zigzag transfer matrix 
(ZTM) (Fig. 10 (dj). The MDEs of the transfer matri¬ 
ces in the infinite-size limit is calculated using the i nfinite 
time-evolving block decimation (iTEBD) algorithirPUMl. 
We denote the bond dimension in the iTEBD calculation 
as Caux- The correlation length Saux of the the statisti¬ 
cal model is plotted against Caux in the upper panel of 
Fig. [TT} When we view the MDE as a ground state 
of a ID chain, we can measure semi-infinite chain en¬ 
tanglement entropy S'aux- Both Saux and 5'aux saturate 
as Caux increases, indicating the uniqueness of the MDE 
and finite-ness of the physical correlation length. Due to 
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Figure 10: (a) The tensor network for the wave function 

of l^'f) is shown. ni, 2,3 are unit vectors along the 30°, the 
150° and the 270° directions, (b) The tensor network for the 
wave-function norm (’Ff |Tf) is shown. When the physical 
degrees of freedom are first contracted, this tensor work can be 
viewed as a partition function of a classical statistical model 
on the auxiliary degrees of freedom, (c) The armchair transfer 
matrix (ATM) on an infinite plane is shown, (d) The zigzag 
transfer matrix (ZTM) on an infinite plane is shown. 



Figure 11: The correlation length Haux of the auxiliary de¬ 
grees of freedom (a) and the semi-infinite chain entanglement 
entropy Saux (b) of the MDE of the ATM are plotted against 
the bond dimension ((aux- As we can see, both quantities have 
saturated at large ^aux. The results for the ZTM are similar. 


the conjugate cycle. We argue that this long range order 
should persist even when the torus is deformed to a finite- 
circumference infinite cylinder, resulting in degeneracy of 
the MDEs of the cylinder transfer matrices. As we in¬ 
crease the circumference L and approach the planar limit, 
the degeneracy should not be lifted due to its topological 
nature. However, in the planar limit, the ATM and ZTM 
each only yields a unique MDE. Therefore, we prove by 
contradiction that [dlf) is not topologically ordered and, 
thus, a spin-1 featureless paramagnet. 


the geometry of the ATM and the ZTM, this result only 
suffices to prove that ) is invariant (without sponta¬ 
neous symmetry breaking) under the translation with a 
doubled unit cell, the SO(3) spin rotation, T, (Tj_ and 
CT||. To show that i? 2 x and fully translation symmetry is 
not spontaneously broken, we calculate the NN spin-spin 
correlation function along the hi , h 2 and h^ (see Figl^. 
Their values agree with each other: {S ■ .S)nn — —0.110. 
Therefore, the full translation symmetry and the 
symmetry are not spontaneously broken. At this point, 
we have confirmed that the wave function I'I'f') preserve, 
without spontaneous symmetry breaking, all the required 
symmetries. 

With the uniqueness of the MDE of the transfer ma¬ 
trices and the symmetry of we now argue that 

the state [Tf) is free of any topological order. This argu¬ 
ment is similar to the one provided in Sec. |III B| We first 
assume that |d>f) is topologically ordered. On a finite- 
size torus geometry compatible with 7?^, the invariance 
of l^'f) under implies that it is an eigenstate of 
ST, where S and T are the topological S matrix and T 
matrix. An eigenstate of ST cannot be an MES along 
either non-trivial cycles of the torus. Consequently, the 
Wilson loops around one non-contractible cycle of the 
torus develops a long range order along the direction of 


C. Spin-1/2 Featureless Paramagnets on the 
Honeycomb Lattice 

In this subsection, we present the slave-spin PEPS 
study for the spin-1 /2 featureless paramagnet on the hon¬ 
eycomb. This subsection contains two parts. In the first 
part, we focus on the slave-spin PEPS construction with 
bond dimension y = 2. Due to the restriction of x, we 
only identify two ’’almost featureless paramagnets”. T 
and at least one of cr_L and a\\ is broken. Numerical stud¬ 
ies verify their symmetry properties. In the second part, 
we propose a slave-spin PEPS construction for a spin-1/2 
featureless paramagnet with bond dimension x = 6. 

1. Spin-1/2 slave-spin PEPS with x = 2 

We start with 3 auxiliary spin-1/2 degrees of freedom 
on each site with each spin-1/2 degree of freedom asso¬ 
ciated to a different link that connects to the same site. 
Then we project the two spin-1/2 degrees of freedom as¬ 
sociated to the same link into a singlet pair. Next, we 
project the 3 auxiliary spin-1/2 degrees of freedom on 
each site into a total-spin-1/2 channel to obtain a physi¬ 
cal spin-1/2 wave function. There are 2 linearly indepen¬ 
dent total-spin-1/2 channels among 3 spin l/2s, resulting 









Figure 12: This figure is the graphical representation of 
Eq[TT] The directly connected lines represent 5-functions be¬ 
tween the indices on their end points, while the lines decorated 
by an arrow represent the singlet projector S defined in Eq. 

The physical 1/2 degrees of freedom are represented by the 
red lines and auxiliary ones by the black lines. 


in the 2 choices of tensors on the site: 

Km = (H) 

Each of the indices fi, a, /3 and 7 run from 1 to 2 rep¬ 
resenting the 2 states of a spin- 1/2 degree of freedom. 
Under the cyclic permutation of the subscripts, they sat¬ 
isfy 


pM _ pZ-n-i/a pM 
a/S'y ^ 




( 12 ) 


Under T and “mirror reflection”, we have 


_ 

jSa.'y 



_ zt-'M 
'^a/ 37 ’ 

_ p'M 

^ a/S-y' 


(13) 


To construct the slave-spin PEPS on the honeycomb lat¬ 
tice, we can either use the same tensor, say , on ev¬ 
ery site or use two different tensors for the 2 sub-lattice 
of the honeycomb lattice. We denote the PEPS associ- 



global spin rotation and . Regarding T, cr± and CTy, 
preserves the combination of any two. ) 

is only invariant under Tcry and crj_. Therefore, we refer 
to them as almost featureless paramagnets. 

Numerical study for - Here, we show numeri¬ 
cally that being an almost featureless paramag- 


net, does not suffer from spontaneous symmetry break¬ 
ing and topological order. With y = 2, we are able 
to perform diagonalization of the ATM and ZTM on fi¬ 
nite circumference cylinder. For the circumference L = 
2,4, 6 , 8 ,10,12, both the ATM and the ZTM has a unique 
MDE in each case. The physical correlation length ^ is 
plotted against 1/L in left panel of Fig. 14 ^ saturates 
to a finite value as 1/L —>■ 0 . Since the transfer matrices 



Figure 13: Different correlation functions including spin- 
spin correlation functions and the triple spin chiralities are 
calculated to verify the unbroken symmetries. Here, ni, 2,3 
are the unit vectors along the corresponding directions. 


are not compatible with certain symmetries, this result 
only suffices to rule out the spontaneous breaking of the 
translation with a doubled unit cell and the global spin 
rotation symmetry. To verify the full translation and the 
All symmetry, we calculated the 3 NN spin-spin correla¬ 
tion functions along three different directions hi, and 
fi-a. Their value agree with each other: {S-S)f^N ~ —1.08. 
Therefore, the full translation symmetry and the are 
confirmed. For the symmetries T(T_l, Tcry and (Tycr_L, we 
calculated the symmetry related triple spi n chiralities: 

for the refer- 


13 


Cijk = {S{fi) ■ (S'(7) X S'(ffc)))(see Fig. _ 

ence the subscripts) and find the following agreement: 


C56I ^716 -C712 C56I, 

Csei = Cyie = -C 712 = -0.05970 ± 0.00002. (14) 


With this agreement plus other checks above, we verify 
that indeed preserves (without spantenous sym¬ 

metry breaking) the following symmetry: the global spin 
rotation symmetry, the translation symmetry, the 


symmetry, tTy(Tj_, Tcr_L and Tcry. 

With the uniqueness of the MDE of the transfer matri¬ 
ces and the invariance under , we can argue, follow¬ 
ing similar strategy in Sec. |IVB and in Sec. IIIB that 
\^K 2 ) topological order. Nevertheless, we also 

calculate the entanglement entropy of the semi-infinite 
cylinder with circumference L (see the right panel of Fig. 


141 from which we obtain an “area law” and the TEE 
0. Therefore, we confirm again the absence of topo- 
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logical order in summary, is indeed an 


almost featureless spin- 1/2 paramagnet. 

Numerical study for \^K 2 ) ~ H^re we present the nu¬ 
merical study for cylinder with circum¬ 

ference L = 2,4,6,8,10,12. Exact diagonalization of 
the ATM and ZTM yields a unique MDE in each case. 
The physical correlation length ^ is plotted against \/L 
in the left panel of Fig. ^ saturates to a finite 

value as 1/L —>■ 0. With these results, we establish the 
absence of the spontaneous breaking of the translation 
with a doubled unit cell, the global spin rotation sym¬ 
metry and (T_l symmetry. We also calculated the 3 NN 
spin-spin correlation functions along three different di¬ 
rections hi, 772 , and 77 . 3 . Their value agree with each 
other: (T-T)nn — —0.076. Therefore, the full translation 
symmetry and the symmetry are not spontaneously 














9 



is mentioned above, both and ) have AF NN 

spin-spin correlation functions. The NNN spin correla¬ 
tion functions are also anti-ferromagnetic: (*S'i*S'5)nnn — 
—0.035 and (<^ 155 )^^^ — —0.041. Therefore, these states 
would be favored by Hamiltonians with AF Ji and J 2 
Heisenberg couplings. 


Figure 14: Left panel: The physical correlation lengths ^ (in 
the unit of lattice spacing) obtained from the ATM (black) 
and the ZTM (red) are plotted against l/L. Right panel: The 
semi-infinite cylinder entanglement entropy is plotted against 
the circumference L. The dash line is a linear fit that extracts 
the “area law” and the TEE 7 . 


2. Spin-1/2 slave-spin PEPS with y = 6 

In the previous subsection, we study the slave-spin 
PEPSs with the bond dimension y = 2 and find two 
almost featureless paramagnets. In this section, we will 
propose a slave-spin PEPS with x = 6 for a fully feature¬ 
less spin-1/2 paramagnet for the honeycomb lattice. We 
start with an auxiliary spin-1 and an auxiliary spin-3/2 
degrees of freedom on each site of the honeycomb lat- 


•- -• ATM 





tice. We put the spin-1 degrees of freedom on 

•- « ZTM 

0 0.1 

0.2 0.3 0.4 0,5 0 

6 



in the wave functions 4>f) studied in Sec. IV B 


l/L 


Circumference L 


Figure 15: Left panel: The physical correlation length ^ (in 
the unit of lattice spacing) obtained from the ATM (black) 
and the ZTM (red) are plotted against 1/L. Right panel: The 
semi-infinite cylinder entanglement entropy is plotted against 
the circumference L. The dash line is a linear fit that extracts 
the “area law” and the TEE 7 . 


broken. Regarding the T<7\\ symmetry, we calculate the 
symmetry related triple spin chiralities: 


using the the spin-3/2 degrees of freedom on each site. 
Both of l^kf) and the AKLT state are featureless. The 
physical spin-1/2 wave function (denoted as |'I'i/ 2 )) is 
obtained by projecting the spin-1 and the spin-3/2 de¬ 
grees of freedom on the same site into a total-spin-1/2 
channel. Since |'Ff) and the spin-3/2 AKLT are both 
PEPS with bond dimension 3 and 2, 1 ^ 1 / 2 ) is a PEPS 
with X = 6. \^i/ 2 ) is, by construction, invariant under 
the global SO{i) spin rotation, the translation, the 
T, cr_L and (T||. 


T crii 

C'rie - > —C'jYi-, 

Crie = - C^x2 = -0.0493 ± 0.0003. 


(15) 


With this agreement plus other checks above, we ver¬ 
ify that ) indeed preserves (without spontaneous 

symmetry breaking) the following symmetry: the global 
SO{3) spin rotation symmetry, the translation symme¬ 
try, the R 2 ^ symmetry, crj_, and Tcry. 

Similarly, the uniqueness of the MDE of the transfer 
matrices and the invariance under implies the ab¬ 
sence of topological order in |'kf/^). Nevertheless, we 
calculate the entanglement entropy of the semi-infinite 
cylinder with finite circumference L (see the right panel 
of Fig. 15) from which we obtain an “area law” and the 
TEE 7 ~ 0. Therefore, we confirm again the absence 
of topological order in ). In summary, 1'!'/^^ ) is 

indeed an almost featureless spin- 1/2 paramagnet. 

To sum up this subsection, the two slave-spin PEPSs 
and ) are both almost featureless param¬ 

agnets on the honeycomb lattice with spin 1/2 per site 
(with different symmetries). Before we end this sec¬ 
tion, we briefly comment on the possible Hamiltonians 
for these two states or for the phases they represent. As 


We perform the iTEBD study for |'I'i/ 2 )- The result 
is shown in Fig. As we increase the bond dimension 
Caux in the simulation, both the correlation length Saux 
for the auxiliary degrees of freedom and the entangle¬ 
ment entropy ^aux of the effective ID chain given by the 
MDE grows linearly. With this result, we cannot distin¬ 
guish the two possibilities: ( 1 ) the transfer matrix has 
a gapless spectrum or ( 2 ) the spectrum is gapped but 
the MDE need to be captured by C > 1200. This leaves 
an interesting question open: whether or not |'I'i/ 2 ) is 
a featureless paramagnet. Note that Ref. [9] proposed 
a featureless wave function I'Feg.) of spin- 1/2 fermions 
at site Filing 1 on the honeycomb lattice. Naively one 
could Gutzwiller project this wave function to produce 
a featureless spin- 1/2 wave function, but they find the 
Gutzwiller projection actually annihilates the wave func¬ 
tion: a spin-1/2 paramagnet is surprisingly evasive. In 
fact, the fermionic wave function I'Feo-) is not even adi- 
abatically connected to our proposal | 4 'i/ 2 ) while pre¬ 
serving the lattice symmetries. That is because the two 
wave functions have opposite eigenvalues under a± on an 
Lx L torus with L odd. Hence, \'^ecr) and | 4 >i/ 2 ) belong 
to different classes of fragile Mott insulatorP^ that are 
separated from each other by a phase transition man¬ 
dated by the lattice symmetries. 
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(a) 


(b) 



Bond Dimension 



Figure 16: The correlation length Haux of the auxiliary 

degrees of freedom (a) and the semi-infinite chain entangle¬ 
ment entropy Saux (b) of the MDE of the ATM are plotted 
against the bond dimension (^aux- As we can see, both quan¬ 
tities haven’t saturated even at the largest bond dimension 
Caux = 1200. 


V. CONCLUSION AND DISCUSSION 

In this paper, we study the problem of construct¬ 
ing featureless paramagnets on the spin-1 square lattice 
and the spin-1 or spin-1/2 honeycomb lattice in 2-1-lD. 
The featureless paramagnets we focus on are short-range 
gapped states with certain symmetries including the 
global 5'0(3) spin symmetry, the translation symmetry, 
the lattice rotation symmetry, the time reversal symme¬ 
try and other possible point group symmetries. Even 
though these featureless paramagnets are not forbidden 
by the LSM theorem or its generalizations, there are 
strong field theoretic obstacles to construct them. Nev¬ 
ertheless, to prove their existence, we use the slave-spin 
PEPS construction to write down ideal wave functions 
for the featureless paramagnets on the square lattice and 
honeycomb lattice. For the spin-1 slave-spin PEPS wave 
functions on the square lattice and the honeycomb lat¬ 
tice, we provide both analytical and numerical studies to 
verify their featureless-ness. For the honeycomb lattice 
with spin 1/2 per site, we present two types of construc¬ 
tion of slave-spin PEPS with bond dimension x = 2 and 
X = 6. With X = 2, we propose and verify two dif¬ 
ferent almost featureless paramagnets with broken time 
reversal symmetry and the mirror reflection symmetries. 
With X = 6, we propose a completely featureless slave- 


spin PEPS wave function. However, due to the lack to 
numerical evidence, the status of the spin-1/2 honeycomb 
model remains ambiguous. 

All these featureless paramagnets are given in PEPS 
forms. Parent Hamiltonians can be constructed such that 
they are the exact ground stated. Sometimes the parent 
Hamiltonians take very complicated forms. It would be 
interesting to find simple Hamiltonians such that their 
ground states are well-approximated by or belong the 
same phase as these ideal featureless paramagnetic wave 
functions. We’ve provided hints to what terms to include 
in the Hamiltonians so that these ideal wave functions 
are favorable. As we’ve shown, for the the case of the 
spin-1 featureless paramagnet on the square lattice and 
the spin-1/2 almost featureless paramagnets on the hon¬ 
eycomb lattice, the NN and NNN spin-spin correlation 
functions are shown to be anti-ferromagnetic. Therefore, 
the Hamiltonians with AF Ji and J 2 Heisenberg cou¬ 
plings is a good direction to look at. 

As we’ve pointed out, the Berry phase patterrP^ in 
the spin-coherent-state path integral provides obstacles 
to the existence of these featureless paramagnets studied 
in the paper. In fact, the Berry phase pattern are ob¬ 
tained in the vicinity of the Neel ordei^^ and the valence- 
bond-solicP. One interesting question would be whether 
it is possible to have a direct second order phase tran¬ 
sition from these featureless paramagnetic phases to the 
Neel order or the VBS order. More generally, one can 
ask what are the classical orders that is separated from 
the featureless paramagnetic phases by a second order 
transition and how one should describe the transition in 
the field theoretic language. In this paper, we obtain 
the ideal wave functions from the slave-spin PEPS con¬ 
struction. In the vicinity of a second transition, if the 
slaved-spin degrees of freedom becomes dynamical, one 
might need a field theory description that is beyond both 
the Landau paradigm and that is a generalized version 
of the deconfined quantum criticality. 
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